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We consider quantum states under the renormalization-group (RG) transformations introduced by 
Verstraete et al. [Pliys. Rev. Lett. 94, 140601 (2005)] and propose a quantification of entanglement 
under such RG (via the geometric measure of entanglement). We examine the resulting entangle- 
ment under RG for the ground states of "matrix-product-state" (MPS) Hamiltonians constructed 
- - , by Wolf et al. [Phys. Rev. Lett. 97, 110403 (2006)] that possess quantum phase transitions. We 

■ find that near critical points, the ground-state entanglement under RG exhibits singular behavior. 
' The singular behavior under finite steps of RG obeys a scaling hypothesis and reveals the correlation 

^ , length exponent. However, under the infinite steps of RG transformation, the singular behavior is 

. rendered different and is universal only when there is an underlying conformal-field-theory descrip- 

' tion of the critical point. 

^3 : 

1—5 , PACS numbers: 03.65.Ud, 03.67.Mn, 64.70.Tg, ll.lO.Gh 

^ : 
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I. INTRODUCTION 

(— I ! 

O • Since Wilson [l|, renormalization-group has been an important tools for theoretical physics, ranging from high- 
^ ' , ' energy physics to condensed matter [2]. It is related to the coarse-graining procedure of the physical system, and 
< from that the transformation of system parameters in the Hamiltonian is derived. Corresponding terms in the 
' Hamiltonian be can determined to be relevant or irrelevant under the scale transformation. The renormalization- 
^ , group transformation on quantum states was recently introduced by Verstraete et al. Q using the representation 
I ^ i ' of Matrix Product States (MPS) (for a review of MPS, see e.g. Ref. [|). Many important quantum states in 
the quantum information theory emerge naturally from the fixed points of this coarse-graining transformation on 

■ states 

^ \ One important property associated with quantum states is their entanglement content. There have been tremendous 

■ advancement on the understanding of entanglement (both bipartite and multipartite) for the past few decades @. 

■ The notion of entanglement has also been applied to many-body systems and especially systems that possess 
quantum phase transitions 3j niainly via bipartite measures (e.g., between two spins or between one subsystem and 
the remaining) such as concurrence [9| and entanglement entropy [lO| . Important insight has thus been obtained, such 
as the connection to conformal theory near criticality [l^l and to the bipartite entanglement [1, [ll| . Furthermore, in 
hig her dimensions than one, the entanglement entropy has been shown to possess an area law for various systems |12| - 

\ Il4| (up to logarithmic corrections in fermions [l2|). Regarding the coarse-graining process, one can thus raise the 
interesting question how the entanglement behaves under the RG transformation on states [151 . As we shall see below, 
one of the measures that appear to be suitable to the discussions of entanglement under RG is the so-called geometric 
measure of entanglement (or simply geometric entanglement) (T6l. [T7|. This measure of entanglement defined with 
$H ' respect to partitions into blocks of consecutive sites has recently been employed by Oriis jl8j and by Botero and 
5^ ; Reznik 

In this paper, we provide an interpretation of the block geometric entanglement, namely that it is exactly the 
entanglement under the coarse graining of the renormalization group transformation on quantum states Q • We apply 
this block-L entanglement (with L being the number of sites in each block) to two spin models constructed by Wolf 
et al. [l^l that possess quantum phase transitions. We find that near critical points, the ground-state entanglement 
under RG transformation exhibits singular behavior. The singular behavior reveals the correlation length exponent. 
However, under the infinite steps of RG transformation, the singular behavior is rendered different and it has no 
universal form unless the critical point can be described by a conformal field theory. 

Let us begin by discussing matrix product states. It was shown by Vidal [21] that any state can be written in the 
MPS form as follows, 

di d2 dm 

l^) = E E - E Tr(4\U[f ..4™1)1pi,P2,-- - (1) 

P1=1P2 = 1 Pm = l 

as long as the matrices Ap'^^'s have large enough dimensions Dk-i x Dk- For any one-dimensional translationally 
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invariant state, all the A's are identical, and the state can be written as 

d 

,P2, ■ ■ ■ ,Pm), (2) 

pi,P2,...,Pm = l 

where A's are D x D matrices with dimension D < d™!"^ . Let us now define the geometric measure of entanglement, 
by considering a multipartite system comprising m parts, each of which can have a distinct Hilbert space. We compare 
this general m-partite entangled pure state jV') to the set of general product pure states, 

\^) ^ ® 10^*^), (3) 

and define the maximal overlap of jV') with the closest product states as follows, 

Ai(V') =max|(</)|V')|. (4) 

The maximal overlap Ai(?/;) reveals the entanglement content of the state 1-0), the larger Ai (?/;), the lower the entan- 
glement of IV')- A quantitative way to define the entanglement content is via 

ii;i(V') = -iogAi(v^)2, (5) 

where the subscript 1 indicates that the product state is composed of product of states of single sites. Note that the 
norm square of the translation invariant matrix product state (with m sites) can be expressed in terms of the operator 
E defined below in Eq. p^ . 

(VIV)-Tr(^'"), (6) 
which, for convenience, may not usually be normalized to be unity, and thus one has to supply this factor in A('0), 

A^W= max MV^. (7) 

We remark that by appropriate partitioning of |0) into various product forms, a hierarchy of entanglement can be 
obtained [13, [H-flll ■ The most relevant kind of partitioning regarding RG is to divide m sites into blocks of several 
consecutive neighboring sites, e.g., L consecutive sites in one-dimension. This leads to what we shall refer to as the 
block-L entanglement, 

£;i(V') = -logAi(^A)^ (8) 

where 

Al(V') = max|($i|V')|, (9) 
with being product states of the block form: 

= , (10) 

and we have implicitly assumed that the total number of spins to (sometimes denoted by N') is an multiple of L. 



II. RG ON QUANTUM STATES 

Verstraete et al. Q considered a quantum coarse-graining procedure by merging two neighboring spins to one new 
block spin, which is described in terms of matrices A's: 

D 

i(f (11) 
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From this a more convenient representation is to choose the new matrix as 

^ a; = \iv\ (12) 

where is the right unitary matrix in the singular- value decomposition 

min(d^,_D^) 

E UI^'^'^K.)- (13) 
1=1 

This keeps the dimension of the Hilbert space of the block to be of size bounded above by . They further introduced 
a more convenient representation of the RG transformation by defining (which we shall thereafter refer to as the RG 
operator or the transfer-matrix operator) 

d d D 

^ = ^Ap®^; = ^ ^ [A^U®[A;]^,\a^i){pvl (w) 

p=l p=l a,p,p,,ii=l 

which is invariant under any local unitary Aq — >■ '^pU!^Ap. Hence, the operator E"^ is invariant under any unitary 
within the block, 

= ^ i(P9)* = ^ A;F' (g) A;F'* . (15) 

pq I 

The RG transformation on the state can then be described by the mapping 

E^E' = E^. (16) 

The above discussions assume that the states being considered are translation invariant. A straightforward extension 
of the renormalization-group transformations on generic quantum states can be made by including the site-dependence 
of the matrices Ap' and E^^K Then one can merge two spins at sites 2fc — 1 and 2k and form a new spin at the fc-th 
site on the new lattice: 

i^f^)^E[Af-il]„,[42'=l],„ (17) 

and 

^^[2fc-l] .^[2fc]^ (18) 

As above, E' is invariant under any unitary on the original two spins. This means that under one-step RG transfor- 
mation, the state \^) transforms to 

1^) ^ 1^') = U[12] (g) U[M] (g) ... (g> U[2k - 1, 2k] (g> (19) 

and the (2fc— l)-th and 2fc-th sites are merged into a single site. The unitaries are now generally site-dependent. 
Under RG of merging two neighboring sites, 

RG2 : IV) ^ IV''). (20) 



III. ENTANGLEMENT OF STATES UNDER RG 



In a similar spirit, though not identical, to the work on multiscale entanglement renormalization ansatz (MERA) 
by Vidal 25], we discard the short-range description and hence the total entanglement should decrease under this 
coarse graining. As the RG defined by Verstraete et al. Q transforms the original state to the state IV'') up to local 
(treating sites 2fc— 1 and 2k as local) unitary transformations [see Eq. ([19])], the natural definition of the entanglement 
after the RG is to determine the one that is minimum among the local equivalence class. In terms of Vidal's MERA, 
the untiaries U[2k — 1, 2k] act as disentanglers that aim to reduce the entanglement between sites 2k — 1 and 2k. The 
merging of the two sites here under the RG defined by Verstraete et al. [3] is done with the same pairs of sites, in 
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contrast to MERA, where the merging is done with, e.g., sites 2k and 2fc + 1 via isometries [25|. In conforming the 
former picture of RG on states, the entanglement after one-step of RG should be defined as follows: 

EG{W})^rmnEGm, (21) 

where the unitary U is of the form U\12\ (g) J7[34] (g) • • • (g) [/[2fc — 1, 2k] (g • • • . The maximal overlap in obtaining the 
entanglement is 



max 



f g) (<?!)W|W|V) =max f ^^(C/^[2fc-l,2fc](/)[2'=-i] ,g,/,[2fc]|\ = max (22) 

Vi=l / 4,'s,U \k=l ' * 



where |$) = |(^[^^l)(g|^[^''l)(g- • •(g|(^[^''~^'^''])(g- • • with |(/)[2*:-i.2fc]^ being any arbitrary state of the two sites 2fc-l and 
2k. This is exactly the block-2 entanglement defined previously in Eq. ([8|) with L = 2. One can continue the merging 
procedure, and arrive at the successive entanglement under RG being equal to El{-4>) with L = 4, 8, . . . , 2', etc. It it 
straightforward to see that for > £^2'' (V") ^ ^2' ("0)1 the total entanglement under RG cannot increase. As 
another irnportant ingredient in the RG is the rescaling of the lattice spacing, we also introduce entanglement per 
block [l^, O, [l9| to reflect this rescaling of length scale (equivalently, system size) in RG, 

£L{i^) ^ (23) 

where N is the total number of sites (usually considered in the limit of large number of blocks, n = N/L — > oo). 
Therefore, E]^ is the total entanglement for the system with merging L sites into one, and correspondingly, is the 
entanglement per block. 

A conclusion from the above discussions is that the entanglement per block (of size 2') is equal to the entanglement 
per site of the Z-time RG transformed state of 1-0), as the RG transformed state is determined up to a 2'-local sites in 
view of the original sites. That is 

£:i(0)=fi(RGf'°^^^(V')). (24) 

This gives a physical meaning to the block-L geometric measure of entanglement. 

Orus [l^l has recently shown that the geometric measure of entanglement defined relative to blocks (of size L) of 
spins can be evaluated via 

MmaxP = max|(fg)f*)^i?^(f (gr*)|, with|r1 = l, (25) 

f 

and 

f = - Km 1 log J^^^ = - logdd^axH, (26) 
n-i-oo n Tr(i?"-^) 

where the important assumption is that the closest product state can be taken to be a product of identical local 
states. For the ground states of the transverse-field XY spin chains, this ansatz has been verified numerically [iGj . 
Let us briefly describe the proof. Consider the ansatz product state to be j^) = j^)®™ and the translation-invariant 
state 1-0) is expressed in the MPS form ([2]). Then the overlap between the two states is 

d 

($10) = Tr(Ap,Ap,...Ap,J(0|pi)(0b2)---(0b™) (27) 

pi,P2,...,Pm = l 

= Tr(S"), (28) 
where the matrix B is defined as 

B^Y.Mm- (29) 
p 

Suppose the largest eigenvalue of matrix B is not degenerate, then 

hm Tr(B™) = A'^". (30) 



5 



Suppose the corresponding eigenvector is r, then the goal is to find 

|Ai| = max|7^i3r| 

r 

= max| ^^r;[Ap]a/3((/>b)r/3|. 

P a,/3 

The maximal overlap becomes (in the limit m — > oo and taking into account normalization 

K{^) = maxmaxl ^^<[AJ„^(0|p)r^ "V|Tr(^")|i/2. 

P a, 13 



Under the assumption that the two maximizations can be interchanged, we get 

A(^) = maxmax^^<[AJ„^(0|p)r^ /|Tr(i?™)|i 

r <p I 

P a, 13 



/2 



The maximization over cj) can be achieved when (using the Cauchy-Schwartz inequality), up to a normalization, 



(31) 
(32) 

(33) 
(34) 
(35) 



This leads to 



Aii^r = max 1 5] r:r^'[ApU[Ap]:,^,r^r;, /|Tr(£;'")| 

p a,a',f},fi' 

= max I Yl <r^'Ec.,c.';P,P'rpr}, /|Tr(£;™)| 
= max (rf(8)f*^£;r(^f*) /|Tr(i;™)|, 



(36) 
(37) 
(38) 



where E is defined in Eq. (|14p as the basis operator for state renormalization. Therefore, the entanglement per site 
becomes 



£i{tp) = lim logA(-i/')^ = lim log 

m— foo 777, m— foo 777, 



(rf (8)f*^£;r'® r*) /|Tr(£;™)| 



(39) 



In the case of entanglement per block of size L, we can repeat the previous derivation by replacing E by E^ and we 
arrive at the expression 



£l{iP) — lirn log 

n— >oo TL 



max 



(ff (g) 7=-*^^;^ f f *) /|Tr(£;"-^)| 



(40) 



Therefore, the entanglement at large block size L — >■ oo depends on the fixed-point property of the operator E. 

IV. FIXED POINTS 

Verstraete et al. Q have also defined the fixed point under the RG transformation on states via 

E°" = lim &. 

1^00 



(41) 



We can use their results on the classification of the fixed point and investigate the behavior of entanglement for generic 
states. They concluded that in the generic case, the largest eigenvalue of E is nondegenerate and both its left and 
right eigenvectors have a maximal Schmidt rank. 



E^ = \<Pi,){<S>lI 



(42) 



6 



and one can always choose 



D 



(43) 



where A, > 0. Then 



{r\{r*\E-'\r)\r*)=J2\,\ 



(44) 



and its maximum is Ai, the largest of {A^}. Furthermore, Tr[(ii^°°)"] — (^X^i Ajj • Therefore, the entanglement per 
block is 



Coo = - hm - log ■ 



Y.tLi A. 



log 



(EvAi 

i=l 



logAi < logD, 



(45) 



where we have defined the normalized Schmidt coefficients Xi = Xi/ A^. This means that if a many-body state can 
be represented by a MPS with dimension D, then the largest entanglement per block at the fixed point is bounded 
above by log!?. We remark that according to Verstraete et al. the entropy of a block of spins is exactly twice the 
entropy of entanglement of i-c, S = —2 Ai log A.; < 2 log D. In short, we have the relation £oo < S/2 < log D. 



V. EXAMPLES 

In this section, we shall warm up by several example states. 
Example 1. AKLT state 



(46) 



Orus has performed a detailed analysis for the AKLT state j23]- The operator E is calculated to be (the matrix 
elements conveniently expressed in the "ket" and "bra" notations, 

E = 3^(100) + |11))^((00| + (111) - 1 • -i=(|00) - |11))^((00| ~ (111) ~ 1 . |01)(01| - 1 • |10)(10|. (47) 



Let \r) = ro|0) + ri|l) and |r*) = rglO) + r^ll), we have 



(r|(r*|^^|r)|r*) = ^ + ^l^. 



Moreover, for a total of n blocks. 



Tr(i;"^) = 3"^ +3(-l 



(48) 



(49) 



So 



Example 2. GHZ state. 



[31 + MOll" 
= - J^-^ 3»U3(^1)»^ = 2 - 



3^ 



(50) 



= 2|00)(00| + 2|11)(11| 



(51) 
(52) 



Then, 



(r|(r*|£;^|r)|r*)=2^(|ro|^ + |rin<2^ 



(53) 
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when either |ro| = 1 or |ri| = 1. Furthermore, 

Tr(i?"^) = 2-2"-^. (54) 

Therefore, 

El ^-\ogj-^=\og2, £l=0. (55) 

Note that in this case, there is a degeneracy in eigenvalues of E, and there seems to be no problem in carrying out 
the procedure. Furthermore, no matter how large L is, GHZ state looks identical (after rescaling the system size), 
and it possesses a total entanglement log 2, and hence vanishing entanglement per block in the limit of large number 
of blocks. 

Example 3. Cluster state. 

Consider the linear-cluster state described by 

The E operator is calculated to be 

^ = 2|00)(-- | + 2|11)(++|, ^2 = 4|$+)($+|, (57) 
where |$+) = (|00) + |ll))/y2. For L > I, 

{r\{r*\E^\r)\r*)^2^-^. (58) 

Furthermore, we have 

Tr(i:"-^) = 2"-^, (59) 

and thus arrive at 

EL=n\og2, £L^\og2. (60) 

We note that the above results hold for L — even; see Example 4 and below. In fact for L = 1, it is known that the 
cluster state possesses Ei = YN/2\ log2, where N is total number of spins 2^. This means that £i — (l/2)log2, 
which is half of f L=ovcn- 
Example 4 ■ Anti- ferromagnetic GHZ state. 



1^) = -(101010101...) + 110101010...)). (61) 

One caution is that the closest separable state for block-1 is either jOlOlOlOl...) or jlOlOlOlO...), neither of them 
being translation invariant. But block- 2 is: e.g., |(01)(01)(01)(01)...). Therefore, there is an even-odd difference. Of 
course, the entanglement is the same as the ferromagnetic GHZ state. To avoid this even-odd effect, we will mainly 
consider the block size L to be even. In the following, we will provide numerical evidence in justifying the ansatz 
we use to calculate the entanglement. We note that for permutation invariant pure states, the ansatz of the product 
states being a tensor product of identical single-site states is well justified; see Ref. [29j . 



VI. ENTANGLEMENT IN QUANTUM PHASE TRANSITIONS WITH MATRIX PRODUCT STATES 

Wolf et al. have recently used matrix product states to engineer quantum phase transitions (QPT) with properties 
differing from standard paradigm ^2^ (e.g., analytic ground state energy and finite entanglement entropy for an 
infinite half-chain), but still with diverging correlation length and vanishing energy gap. Since the ground state 
depends on the system parameter (g) continuously, we shall investigate the ground-state entanglement properties of 
the corresponding models Wolf et al. considered and determine whether the ground-state entanglement can be a 
telltale of the corresponding critical points. 
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FIG. 1: Entanglement density £L{g) of the ground state of the Hamiltonian (|63p for (a) A*' = 10 spins and L = 1 (top panel); 
labels "a", "b", and "c" denote the results from the identical, alternative, and arbitrary ansatz states, respectively; and (b) 
L = 2 (bottom panel); here labels "a" and "b" denote the results from the identical and arbitrary ansatz states, respectively; 



A. A spin- 1/2 model 

The state that is represented by 
is the ground state of the following Hamiltonian [20j] 

H = j2 2(5' - 1)'^- ^r+i - (1 + + {9- cTr+i<+2. (63) 

i 

At g = 0, it is a GHZ state, whereas at g = —1 it is a cluster state. 

We first check to what extent the ansatz states for deriving the formula for £l can be justified. In Fig. [5^, 
we compare the numerical values of the supposed entanglement density (£1) with iV = 10 spins using three different 
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-1 . 5 



FIG. 2: Entanglement behavior of the ground state of the Hamiltonian (|63|) . Upper panel: Entanglement £oo{g)- Lower panel: 
Entanglement derivative d£oo{g)/dg. 



product states (identical, alternating, arbitrary) 

l<fa) = |</>^^), (64) 
\^b) = <E) \(l)2) \<Pl) \<p2) ® ■ ■ ■ , (65) 

l<i>c) = (^l^\4>^)■ (66) 

We see that the results of using the ansatz (identical) to calculate £i are correct only for g > —0.5. However, 
the results of using (alternating) are always as good as those of using j^c) (arbitrary). This suggests that for 
L = even, an product of blocks of even number of spins is a good ansatz. This is indeed the case for L = 2. In 
Fig. [1}d, we compare numerical results for £2 with the same number A'^ = 10 of spins between two different ansatz 
states (identical and arbitrary) 

l<i>a) = 10^/^), (67) 
|*b) = 1012)® |'/>34) <»•••«) |0iV-l.Ar). (68) 

It is clear from the plot that they give identical results, thus supporting the use of product of identical single block 
states (with even block size L). 

Next, we calculate the analytic expression for the entanglement per block, focusing mainly on L = even. We begin 
by noting that the operator E is 



2\ 



^=looo o\=a + 9rm{gg\ + 2\n){+ + l (69) 

ail 1 
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Sl 

. 8r 




-3-2-10 1 ^ 

FIG. 3: Entanglement behavior Ehig) of the ground state of the Hamiltonian (|63|) for L = 1,2,4,8 from bottom to top. Note 
that in the top panel for L = 1 only the range g > —0.5 is shown. The dashed curve represents the results of £00(9) ■ 



where \g) = (|0) + g\l))/ \/l + g^. We can evaluate the n-th power of i?, 



/GO 


{l-gY 



\0 








(70) 



where 



U = 



(9 -g g\ 
-l-.g -1 


V 1 



1 
1 1/ 








1 

1+9 


1 

1+9 












1 







1 








1 






2 


I 




1 


1 




2ff 


1+9 


1+9 





(71) 



We can then calculate £l{9)- In this case, we shall see that even the behavior of the entanglement £ao{g) will exhibit 
singularity across the critical point 5 = 0. We calculate the entanglement £oo{g) to be 



r log 2 + log(l + g)-2 log (1 + for g > 0, 

foo(ff) = I 0, for 5 = 0, 

[ log 2, for g < 0. 

The derivative of £oo{g) is discontinuous across g = and exhibits divergence as <? — 0+, 

^^00(5) 



(72) 



dg 



= -g-''\ g^0+. 



(73) 



Therefore, the critical point is reflected by the property that the fixed-point entanglement £00 has divergence behavior 
in its derivative with respect to g. 
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FIG. 4: Nearest-neighbor Concurrence behavior C{g) of the ground state of the Hamihonian (|63|l . The use of concurrence to 
infer critical points may incorrectly identify (7 = 1 as a critical point. 



We can also study the entanglement £l for finite L and learn how the entanglement varies under state-RG trans- 
formation. The entanglement for L finite turns out to exhibit more feature (see Fig. 



log2 + log(l + g) - 2 log (1 + Vff), for .9 > 0, (1 l-.gr > (1+5) 
1 + (1-5)^(1 +.9)-^+. 9(1 +.9)^-2|l 



log 2 - log 
log 2 - log 



^(1 



|g| 



(I + I.9I) 



forg>0,(l+g)^|l-5 
for g <0. 



V5 < (1+5), (74) 



The singularity at the critical point 9 = is obvious; see Fig. [3] In particular, we find that 



dg 



L 



1 dSLig) 



dg 



-{L- 



8' 



(75) 



As g decreases from a large value, the entanglement decreases to zero at g — 1, which is paramagnetic state with all 
spins pointing in the x direction. It then rises to a local maximum as g further decreases, and afterwards decreases 
to zero at the critical point. This intermediate region becomes smaller as L increases and is washed out at the fixed 
point. The cusp for L = 2 at 9 = —1 reflects a highly entangled state, which turns out to be the one-dimensional 
cluster state. The entanglement slowly decreases as g become more negative. 

As a comparison, we show the behavior of the nearest-neighbor concurrence (soj in Fig. 21 There are singularities 
at g ~ and 5 = 1. The use of concurrence to infer critical points may incorrectly identify 5 = 1 as a critical point. 



B. A spin-1 model 

This example is a spin-one scenario. The state represented by 

Ao = -(7z, Ai = cr+, A^i = ga^ , 
is the ground state of the following spin-1 Hamiltonian 

H2 = Y.^2 + g-')SA+i+2{SA+if + 2{A-g''){S, 

i 

-{g + 2f{StSt^,f -f .9(5 + 1){StSt+,jA+^- 



(76) 



z\2 



(77) 



For g — ±2, the GS is the AKLT state. For g ±00, the GS is the Neel GHZ state. The critical point is at 
g — gc — 0, where there is a diverging correlation length. 

Here, we also check to what extent the ansatz states for deriving the formula (|40| for £l can be justified. In Fig.[5K, 
we compare the numerical values of the supposed entanglement density (fi) with = 4 spins using three different 
product states (identical, alternating, arbitrary) 



l^h) = \4>l) «) 102) «) {(l^l) «) \4'2) (E) ■ 
l$c) = ®ili|0,). 



(78) 
(79) 
(80) 
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FIG. 5: Entanglement density £]^^i(g) of the ground state of the Hamiltonian (|77|) (a) for A*' = 4 (top panel); labels "a", "b", 
and "c" denote the results from the identical, alternative, and arbitrary ansatz states, respectively; (b) (bottom panel) for 
A'' = 4, 6, 8 and the analytic expression from Eq. (|85]) . 



0.8r 




-4 -2 2 4 ^ 

FIG. 6: Entanglement Shig) for the ground state of the Hamiltonian (|77l) for block sizes L — 1, 2, 4, 8 from bottom to top. 
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We see that the results of using the ansatz |$a) (identical) to calculate Ei are correct only for \g\ being small. For 
larger \g\ the ansatz \^a} gives incorrect results. However, the results of using |$;,) (alternating) are always as good 
as those of using |<i>c) (arbitrary). Similar to previous model, this suggests that for L — even, an product of blocks of 
even number of spins is a good ansatz. We note that for finite N there is a region near g — that the entanglement 
behaves quadratically with g. This region, however, shrinks as N becomes larger, as illustrated in Fig. [SId. In the 
thermodynamic limit (N — > oo) this region is expected to disappear and the behavior of entanglement near g becomes 
linear (see below). Furthermore, even though the formula for is expected to work for L = even, the expression 
obtained by taking L = 1 appears to be close to the entanglement density obtained for finite N, except the small 
quadratic region (due to finite system size). 

Now, we derive the analytic expression for ground-state entanglement per block. First, we can construct the 
renormalization operator 



E = 



which can be transformed into the diagonal form 



E^U 



/I 
0-100 

0-10 

Vl 1 






1-.9 




(81) 



(82) 



where 



U = 




From this, we can get the £ooig), 









1 





1 





1 








2. 










2s 



Eooig) 



(83) 



log 2, for 5 7^ 0, \g\ < oo, 
0, for g — 0, \g\ — ?> oo. 



(84) 



The critical point in the limit of the fixed-point L — oo is reflected only by the discontinuity of the entanglement. 
The more interesting analysis comes to the case L is finite. For simplicity, take L to be even. For \g\ ^ 1, £l{9) = 0. 
We find that £l(— ff) = ^Lig), as can be seen from the symmetry of E under g — For finite values of g, we have 



^L(.g) 



log 2 
log 2 



log 
log 



l + (l+g)-^(5-l)^], for2<g 
l + (l + g)-^], for0<5<2. 



We find that there is a discontinuity in the derivative across the critical point. 



dSUg) 



dg 



0+ 



dg 



= -L. 



0- 



(85) 



(86) 



The critical point is revealed by this singularity (i.e., the discontinuity). In addition to the critical point, the AKLT 
state aX g = ±2 is signalled by the cusp in the entanglement for L = 2 and 4. The weaker singular behavior at g = ±2 
suggests that it is not a critical point. 



VII. ENTANGLEMENT NEAR CRITICALITY 



Near quantum critical points, the correlation length generally scales as 

^{g)-\g-gc\-r 



(87) 



For any density functions, such as free energy density and the entanglement density, when there is no logarithmic 
singularity we expect the scaling behaves as (i.e., the scaling hypothesis, see, e.g. Ref. [2|) 

£dg)-{i/L)-\ 
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where d is the dimension of the system. We then expect that (unless d£{g)/dg behaves as ~ log \g — gd) in general 

d£L{g) 



^ ±L%-gd^'''\ (89) 
dg s? 

In the case of one dimension d — 1 we consider here, the discontinuity of the entanglement derivative that we have 
found for both models (with L finite) implies that we have v ~ I, which is indeed the case. Furthermore, the 
discontinuity is proportional to L, which is expected as it is now the smallest length scale. 

In the case of transverse-field XX model, near the critical field gc — ^ the entanglement density £i has been shown 
to behave as 



Si^'-^VT^, forg^l-, (90) 

and the entanglement derivative behaves as 

dS^"" 1 



dg VI -5' 



fovg^l-. (91) 



This is consistent with d — 1 and = 1/2 for the XX model. 

These results need to be compared to the case of one-dimensional transverse-Ising spin chain, where the divergence 
of the entanglement derivative is logarithmic, 

^^^^^^-^log|5-l|, 3^1, (92) 

where g is the ratio between the external field and the spin- spin coupling and 17 = 1 is the critical point. In this case 
finite-size scaling needs to be performed to determine 1^ = 1 |16j . 

We remark that the divergence behavior in £l can be rendered different as L — )■ 00 (i.e., much larger than the 
correlation length ^), so the scaling hypothesis is no longer expected to hold. The two MPS models discussed above 
already illustrate this: £finitcL and £00 have different singular behaviors near critical points. For the transverse-Ising 
model, the results of Onis [18] showed that near g = 1, the entanglement £l with block size L 3> ^ behaves as 

£oo{9) - ^\og^ - ^\og\g - (93) 

which gives (where c is the central charge of the conformal field theory) 

In brief, near criticality the singular behavior of entanglement under finite steps of RG reveals the critical point 
and the associated correlation-length critical exponent v. However, under the infinite steps of RG transformation, 
the singular behavior is rendered different, which seems to be non-universal for critical points not describable by 
a conformal field theory, such as the two MPS models discussed here [23|. But it is universal for critical points 
describable by a conformal field theory such as the transverse-field Ising model. 



VIII. FIDELITY MEASURE 

As a comparison, we compute the logarithmic fidelity [3l| between ground states at two different system parameters 
(.9 = 91 and g ^ ga), 

/(,„,.) . hm -log ,,, lit^yif^^J/^ (95) 

m->oo m {V(gij\v(9i)){v(92j\w(92)) 

where m is the number of sites (or spins). It has been shown that the fidelity measure can reveal the underlying 
critical point. The inner product between two matrix product states 

di d2 d^-fi 

W ^ - J2 TiiAp,Ap,---ApJ\pi,p2,--- ,p.m), (96) 

,P2, ■ ■ ■ ,Pm), (97) 
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FIG. 7: Fidelity measure f{g\,g2) between ground states of the same model at two different system parameters g\ and 32- It 
turns out this both Hamiltonians (|63|l and (|77|l have the same plot. 



is straightforward to evaluate: 



{(j)\^) = TYiW- 



(98) 



where W = X^p^ ^Pk ® "^^^ fidelity measure for both models has been calculated by Cozzini, lonicioiu and 
Zanardi for arbitrary m |32| . Here we compare with the results in the limit of to — cx). Interestingly, for both models, 
we obtain the same expression 



1(91,92) 



(H 


-\gi 


)(! + 


92I) 




1+ 


3192 




(H 


-Isi 


)(! + 


92I) 



for gi52 > 0, 
for gi52 < 0. 



(99) 



The results are shown in Fig. [71 The critical point (gi = 52 = 0) is at the intersecting point of four different regions. 
As seen from the figure, there is singularity across the lines (0, (72) and (gi, 0), and (0, 0) is the intersection of all these 
four lines. 



IX. CONCLUDING REMARKS 



We have considered the entanglement of states under the renormalization-group (RG) transformations and apply to 
the ground state of "matrix-product-state" Hamiltonians constructed by Wolf et al. For these models, the entangle- 
ment entropy of a consecutive L spins does not scale logarithmically with L. Furthermore, the use of concurrence for 
one of the models can lead to a spurious critical point. Using the geometric entanglement under RG, we have found 
that near critical points, the ground-state entanglement exhibits singular behavior. The singular behavior under finite 
steps of RG obeys a scaling hypothesis (similar to free energy) and reveals the correlation length exponent. However, 
under infinite steps of RG transformation, the singular behavior is rendered different. It is universal only when there 
is an underlying conformal-field-theory description of the critical point. Along the way, we have provided an upper 
bound for entanglement per block, which is log I?, where D is the dimension of the matrices in MPS. This also shows 
the more complex the ground state is, the larger the dimension of the representative matrices for MPS we need to 
use. 

We conclude by posing the question whether there is any significance to the following function 

Pi9,L).§M. (100) 
alogL 

At least it is a quantity showing how the entanglement changes under the RG scale transformation; see Figs. [3] and ID 
Under RG (i.e., as L increases) certain singular but non-critical features get washed or smoothed away. But the singular 
behavior near criticality persists for large L. The critical points of the two Hamiltonians ([55)) and (|77l) actually have 
different fates approaching fixed points. The former becomes an isolated point in the entanglement, whereas the latter 
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become an algebraic singularity in the entanglement derivative, albeit both singularities are rendered different from 
finite L. 
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